Supersymmetric contribution to the CP asymmetry of B — > J/ip<fi in the light of recent 

B s — B s measurements 
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We derive new model independent constraints on the supersymmetric extensions of the standard 
model from the new experimental measurements of B s — B s mass difference. We point out that 
supersymmetry can still give a significant contribution to the CP asymmetry of B 3 — * J/ip(j> that 
can be measured at the LHCb experiment. These new constraints on the LL and RR squark mixing 
severely restricted their possible contributions to the CP asymmetries of B — » cj>K and B — > n'K. 
Therefore, SUSY models with dominant LR flavor mixing is the only way to accommodate the 
apparent deviation of CP asymmetries from those expected in the standard model. Finally we 
present an example of SUSY non-minimal flavor model that can accommodate the new AMb s 
results and also induces significant CP asymmetries in B s — + J/4><t>, B — > <j>K and B — > r/'K 
processes. 
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1- Recently, the D0 1] and CDF @] collaborations 
have reported new results for the B s — B s mass difference: 



17 ps~ L 
A Mb 



< AM Bs < 21 ps~ 



90% C.L. (£>0), 



17.33±°;g ± 0.07 ps" 1 {CDF), (1) 



which seems consistent with the Standard Model (SM) 
predictions. In fact, the estimation of the SM value for 
AM Bs contains large hadronic uncertainties. The B® — 
B® mass difference is defined as AM Bs = 2A4i2(B s ) = 
2\(B° s \H£ f B = 2 \B° s )\, where H^ B = 2 is the effective Hamil- 
tonian responsible for the AB = 2 transition. In the 
SM, H^ B=2 is generated by the box diagrams with W- 
cxchange. The best determination for AM| M can be ob- 
tained from a ratio to the AMj^ 1 in which some QCD 
corrections as well as t quark mass dependence are can- 
celled out 



AMff " MB~ d 



b bJ Bs 



V',, 



Br 



(2) 



where M Bd = 5.28 GeV and M Bs = 5.37 GeV and the 
lattice calculations lead to BB s f Bg /(BB d f Bd ) — (1.15 ± 
0.06±oZ) 2 S- Since the B d - B°d oscillation is mostly 
saturated by the SM contributions we can assume 
that AAffM = AAf° x d P = (0.502 ± O.QQ^ps- 1 . Finally, 
| | 2 / 1 V* | 2 can be given as a function of the angle 7 
of the unitary triangle of Cabibbo-Kobayashi-Maskawa 
(CKM) mixing matrix. In Fig. 1, we present the allowed 
range of AM§ M in terms of the angle 7 (measured from 
a pure SM process). Here we assume that |V c t,| and \V u b\ 
are free of new physics and can be determined by the 
SM contribution to the semileptonic decay. Also, it is 
assumed that the angle (3 is given by f3 SM , measured 
from Bd — > J/tpK s . As can be seen from this figure, the 
new bounds on AM Bs impose stringent constraints on 
the values of 7 . The lower bound of D% result excludes 
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FIG. 1: Allowed region of AMb,, in the SM, as function of 
the angle 7. 



values of 7 SM > 70°. It is worth mentioning that the best 
fit for 7 SM and AM% M , according to UTfit group is given 
by |: 



7 



SM 



61.3 ± 4.5, AAf| M = (17.45 ± 0.25)ps~ 



and according to CKMfitter group is given by 



7 



SM 



59. 



IX, AM B M = 17.3 



^0.49 

-0.20- 



(3) 



(4) 



Therefore, it is expected that the experimental mea- 
surements in Eq.Q provide important constraints on 
any new physics beyond the SM Jjj. In this letter, we 
study the constraints imposed on supersymmetric model 
due to these experimental limits. We derive model inde- 
pendent bounds on the relevant SUSY mass insertions. 
Then we analyze the implications of these constraints on 
the supersymmetric contribution to the CP asymmetry 
in B s — > J/V"^ process. Finally, we consider the SUSY 
non- minimal flavor model studied in Ref. 8], as an ex- 
ample for SUSY model, that can accommodate the new 
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AMb s results and also induces significant CP asymme- 
try in B — > J/tpcf) which can be measured at the LHCb 
experiment. 

2— In supersymmetric theories, the effective Hamilto- 
nian H^ B=2 receives new contributions through the box 
diagrams mediated by gluino, chargino, neutralino, and 
charged Higgs. It turns out that gluino exchanges give 
the dominant contributions 9]. The most general ef- 
fective Hamiltonian for AB = 2 processes, induced by 
gluino exchange through AB = 2 box diagrams, can be 
expressed as 



3 

# e A / =2 = E ^(M)Qi(M) + E Ci(n)Qi(») + h.c. , (5) 

»=i 



where Cj(/i), Ci(fj) and Qiijj), Qi{n) are the Wilson co- 
efficients and the local operators normalized at the scale 
mi, respectively, which can be found in Rcf. 9] . As in the 
Bd system, the effect of SUSY can be parameterized by 
a dimcnsionless parameter r s and a phase 26 s defined as 
follows: 



I M 12 (B S ) 



(6) 



(B° S \H C ^\B°) 



M SM 



12 T A4 12 



B s mass difference is given by 
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where Mi 2 {B s ) 
Thus, the total B s 
AM B , = 2\M\f{B s )\ r 
sertion approximation, the gluino contribution to the am- 
plitude of B s oscillation is given in terms of the ratio of 
the gluino mass to the average squark mass, x = m^/m^, 
and the down squark mass insertions between second and 
third generations, [640)23, where A and B stand for left 
(L) or right (R) handed mixing. A general expression for 
R s — M^/Myf has been given in Ref.Q as follows: 



R s = a 1 {m 5 ,x)[{Si L )l 3 + (6i R )l 3 ]+a 2 (m 5 ,x)l{Si R )l 3 
+ (S d R Lf2s\+^(m S ,x)[(6 d LR ) 23 (S d RL h3} 



(7) 



where the coefficients |oi| ~ 0[1), \o, 2 \ < \o>3\ < |c*4| — 
0(100). For instance, with mg — 300 and x = 1, one 
finds 



R s = 7.2[[6 d L ) 2 23 + [6 d RR ) 2 23 ] 
205.7[[6 d R ) 23 [S d RL ) 23 ] 



- U9.8[[6 d R ) 2 3 + (S d RL ) 2 23 ] 
803.8[(S d LL ) 23 (6 d RR ) 23 }. (8) 



Note that r 2 s = | 1 + R s |. From the experimental up- 
per bound on AAIb* in Eq.Q, one can derive an upper 
bound on the mass insertions involved in Eq. ■ In order 
to find conservative upper bounds, we set the SM contri- 
bution to its best fit value, namely AM B M = 17.5 ps -1 . 
In this case, the \R S \ should satisfy the following bound: 



\Rs 



(AMbJexp 



(AM Bb ) 



1 



SM 



< 



4/17. 



(9) 
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\{st L u\ 


l(^«)23| 






0.25 


0.074 


0.035 


0.018 


0.014 


0.5 


0.11 


0.037 


0.024 


0.015 


1 


0.17 


0.04 


0.032 


0.016 


1.5 


0.27 


0.43 


0.039 


0.017 


2 


0.46 


0.046 


0.046 


0.018 



TABLE I: Upper bounds on |(di s )a3|, {A,B} = {L,R} from 
AM Bs < 21 ps -1 for nig = 300 GeV. 



It is worth mentioning that if one assumes that AM B M ~ 
21 ps" 1 , the above bound remains valid. In table 1 we 
present our results for the upper bounds on |(<5^ B )23| 
mass insertions from their individual contributions to 
B s — B s mixing for rriq — 300 GeV and x varies from 
0.25 to 2. As can be seen from Eq.|J7|) that the con- 
straints imposed on the mass insertions are symmetric 
under changing L «-> R. Therefore, we present in table 1 
the upper bounds on one combination of the mass inser- 
tions. 

Three comments on the results of table 1 are in or- 
der: 1) the constraints obtained on \{5f jL( ^ RR ^) 23 \ are 
the strongest known constraints on these mass inser- 
tions, since other processes based on b — > s transition, 
like B — > A s 7, leave them unconstrained 0]. In fact 
with these constraints, one can verify that the LL(RR) 
contributions to B — > cj)K, B — ► rj'K and B — > ttK 
are diminished and become insignificant. Therefore, LR 
contribution remains as the only candidate for saturat- 
ing any deviation from the SM results in the CP asym- 
metries or branching ratios of these processes 0- 2) 
The upper bounds on LR(RL) mass insertions from the 
B s — B s are less stringent than those derived from the 
experimental limits of the branching ratio of B — > X s "f 
|l0|. 3) The combined effect of (<5f L )23 and (S RR ) 23 is 
severely constrained by AMb 3 ■ However, the lowest value 
of (S d L ) 2 3, that can be obtained in the minimal SUSY 
model with universal soft SUSY breaking terms, is of or- 
der A 2 ~ O(10~ 2 ). Therefore, it is clear that models 
with large RR mixing would be disfavored by the AMb 3 
constraints, consistently with the previous conclusions 
reached by using the mercury electric dipole moment 
(EDM) constraints |ll|. Indeed, with a large (S RR ) 2 3 one 
may induce a large imaginary part of the mass insertion 
(S d R ) 22 which overproduces the mercury EDM. This also 
implies strong constraints on the right squark mixings. 

3— The B s — * J/4'4 1 decay is accessible at hadron 
colliders where plenty of B s will be produced. It is, 
therefore, considered as one of the benchmark channels 
to be studied at LHCb experiment. The final state of 
B s — > J /ip4> is not a CP eigenstate, but a superposition 
of CP odd and even states which can, however, be dis- 
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entangled through an angular analysis of their products 
. This angular distribution yields to a tiny direct CP 
violation. Thus, the CP asymmetry of the B s and B s 
meson decay to J/ip'fi is given by 



(*) 



= r(B° a (t) -> j/^0) - r(i? s (t) -> j/^) 
r(S«(t) - W) + r(s°(i) - J/v^) 

= sin(AM Bs i), (10) 



where Sj/^^ is the mixing- induced CP asymmetry. In 
the SM, the mixing CP asymmetry Sj/^^ is given by [9J 

Sj,^ = sin2/3 s SM = sin [2 arg(V tb V t * s )] 

~ -2A 2 r;~O(10- 2 ). (11) 

Such small CP asymmetry in the SM gives the hope that 
if a sizable value of Sj/^^ is found in future experiments 
(in particular at LHCb experiment), then it would be an 
immediate signal for a new physics effect. 

In the presence of SUSY contribution, the CP asym- 
metry Sj/^ is given by 



, J/il>(j> 



sin 2/3' 



cir 



in(2/3 s SM + 20 s ) 



(12) 



where 9 S is given in Eq.JBJ) as 28 s = arg(l + R s ). There- 
fore, the value of Sji^ depends on the magnitude of R s 
which, as emphasized above, is constrained from AMb 3 
to be less than or equal 4/17. In this respect, it is easy 
to show that the maximum value of Sj/^^ that one may 
obtain from SUSY contributions to the B s — B s mixing 
is given by 



0.24. 



(13) 



It is important to note that due to the stringent con- 
straints on (S'[ R )23 from b — > S7: \{S'[ R )23\ < 0.016, the 
LR (RL) supersymmetric contribution to Sj/^^ is very 
restricted. It implies that Sj/^^ < 0.02, which is too 
small to be observed at the Tevatron or the LHC. There- 
fore, the LR and RL contributions can not provide sig- 
nificant contribution to B s mixing or to the mixing CP 
asymmetry of B a — > J/ ip(f>. 

On the other hand, the LL and RR mass insertions 
can generate sizable and measurable values of Sj/^p^. 
For instance, {5 d LL(RR) hz ^ 0.17 e i7r / 4 yields to R ~ 
0.24 e"/ 2 which implies that sin2/3 s ~ 0.24. However, 
as mentioned above, it is important to note that since 
the minimum value of the mass insertion (Sf^ L )23 is of 
order 10 -2 , thus, in case of SUSY models with large 
right-handed squark mixings, i.e., (S RR )23 ~ 0.17, one 

O(10 _1 ) which may ex- 



finds that ^0^)23(^)23 
ceed its upper bound presented in table 1. Therefore, 
in this scenario, contributions from both (1^^)23 an d 
(^ii)23(^flii) 2 3 should be considered simultaneously in 
determining the AMb s and sin2/3 s . 



4— We now consider the impact of the AMb, con- 
straints derived above on the mixed CP asymmetries in 
Bd — > <j)K and Bd — > rj'K processes, which at the quark 
level are also based onii->s transition. The BaBar and 
Belle results for these asymmetries lead to the following 
averages: 



S, 



cbK 



0.47±0.19, 



S V , K =0.48 ±0.09, (14) 



which display about ler and 2.5cr deviation from the SM 
predictions, respectively. 

The SUSY contributions to the decay amplitudes of 
Bd — > <f>K and Bd — > rj'K are given by [13j 



i=l 
12 

^ml^^Kf^HiWXQ-Q), (15) 

* i—1 



,G 



where the C are the corresponding Wilson coefficients to 
the local operators of b — > s transition. C a s functions of 
the mass insertions (S'[ L )23 and (8 LR )23 and C as func- 
tions of {5 RR )23 and (S RL )23 can be found in Ref . [l^| . 
Here the QCD factorization mechanism is adopted to de- 
termine the hadronic matrix elements and as in Ref.[l3| 
they can be parameterized in terms of the parameters 
Hi((f>) and Hi{rj') which are given in Ref. 01 . In terms of 
SUSY contributions, the CP asymmetry Sm^k can be 
written as [l4| 



s <t>(n')K = sin 2(3 + 2\R <f , {v ,- ) \cosS (t>W ) sinO^^ cos 2/3, 



(16) 



where R, 



<t>W) 



/ A su 



SY 



4>(r,')K 

in order to derive S, 



and &<j>(r)<) is the strong phase. 
(p(v')K t° war( i their central 



arg 
Thus 

values of the average experimental results in Ea. (|14f> . 
\R(p(ri') \ ^ 0.2 should be satisfied. For a gluino mass and 
average squark mass of order in — rrig — 500 GeV, one 
finds 

R tj> = - O.Ue-* 01 (5 d LL )23 - l27e- t0m (S d LR )23 + L^ R, 

(17) 



and 



R n , = -0.07e i0 - 24 (6 d LL ) 23 - &4(5 d LR ) 23 - L <-> R. 



(18) 



It is now clear that the AMb 3 constraints play a cru- 
cial role in reducing the LL and RR contributions to 
the S<j,ini\K- By implementing the bounds in table 1, 
one can easily observe that LL(RR) contribution leads 
to iR^i-jl ~ O(10~ 2 ) which yields a negligible effect on 
S^ij and one can safely conclude that the LL and RR 
mass insertions can not provide an explanation to any 



4 



deviation in Sm^-j results. On the other hand, the contri- 
bution of (S'[ R )32 is less constrained by AMb, and large 
effects in \R^( v /)\ that could drive S^^k toward 0.4 can 
be achieved. 

5— The above results show that Sj/^^ and S^i^x are 
dominated by different mass insertions: LL and LR/ RL 
respectively. As emphasized in Ref. 8] , these two mass in- 
sertions can be enhanced simultaneously in SUSY models 
with intermediate/large tan/3 and a simple non- minimal 
flavor structure, where the scalar mass of the first two 
generations is different from the scalar mass of the third 
generation. In particular, let us consider the following 
soft SUSY breaking terms are assumed at the GUT scale. 



phase of /x set to zero to overcome the EDM constraints. 



Mi = M 2 = M 3 = M 1/2 , 
M% = ml, 



Ml 



A u = A d 



-Hi 



A e 



'H-2 



Mi 



(19) 



9 2 
a m - 



The parameter a measures the non-universality of the 
squark masses. It is worth mentioning that the EDM 
constraints on the CP violating phase 4>a of the trilinear 
coupling is less sever than the constraints imposed on the 
other SUSY CP phases and can be of order 0(0.1) 0. 

Using the relevant renormalization group equations, 
one can explore these parameters from GUT scale to 
the electroweak scale, where we impose the electroweak 
symmetry breaking conditions and calculate the squark 
mass matrices. Then we determine the numerical values 
of the corresponding mass insertions. For instance, for 



5, tan/3 = 15 and m~„ 



500 GeV, one finds 



that |(5f z )23| - 0.18 which leads to AM Bs ~ 19 ps" 1 . 
Also with a proper choice for the phase 4>a, one can get 
arg[(<5L L )23] — 0.7 which implies that Sj/^,^ ~ 0.1 which 
can be measured by the LHCb experiment. Note that in 
this scenario the phases of the mass insertions are due to 
a combined effect of the SM phase in the CKM mixing 
matrix and the SUSY CP phase 4>a- However, for the LL 
mass insertion the main effect is due the CKM phase, see 
Ref.Q for more details. 

Concerning the mass insertion (S d R )2 3 , it is expected 
to be negligible due to the universality of the trilinear 
couplings. However, with intermediate/large tan/3, the 
double mass insertion is quite important and it gives the 
dominant effect as follows 



(20) 



rrib(Ab —fi tan /3) 



where (S d R ) 33 l= — 
blc > ( (5 Li?)23 off is given by 



Since (8^)23 is negligi- 



(8 d LR ) 23s{! ~ (8 d LL ) 23 -± tan/3. 

m 



(21) 



The parameter /x is determined by the electroweak condi- 
tions and it is found to be of order the squark mass. The 



Since (8 d LL . 



23 



0.18, the value of {8j jR ) 23ett is of order 



10~ 2 which is sufficient to reduce the CP asymmetries 
S^k an d S„'k from the SM result sin 2/3 ~ 0.7 to their 
central values of average experimental results. 

6— To conclude, We have considered the supersym- 
metric contributions to the B s — B s mixing. We derived 
new model independent constraints on the magnitude of 
the mass insertions (8 AB ) 23 , where {A, B} = {L,R}, 
from the new experimental measurements of AM Rs ■ We 
showed that by implementing these constraint, the SUSY 
contribution, through the LL mixing, can enhance the CP 
asymmetry of B s — > J/ip^) up to 0.24, which can be ob- 
served at the LHCb experiment. We also emphasized that 
the new constraints exclude the SUSY models with large 
RR flavor mixing and severely restrict the LL contribu- 
tions to the CP asymmetries of B — > <f)K and B — > rfK. 
Therefore, SUSY models with dominant LR flavor mixing 
is the only way to accommodate the apparent deviation 
of CP asymmetries from those expected in the standard 
model. Finally we studied an example of SUSY non- 
minimal flavor model and intermediate/large tan/3. We 
showed that in this model the new AMg a results and 
also the CP asymmetries in B s — > J/ip(f>, B — > cj)K and 
B — > rj'K processes can be simultaneously saturated. 
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